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We study the problem of critical slowing-down for gauge-fixing algorithms (Landau gauge) in SU(2) lattice 
gauge theory on 2 and 4 dimensional lattices, both numerically and analytically. We consider five such algorithms, 
and we measure four different observables. A detailed discussion and analysis of the tuning of these algorithms is 
also presented. 



1. INTRODUCTION 

Here we report on the status of our study |Q,D 
of Landau gauge- fixing algorithms. The efficiency 
of these algorithms is of key importance when 
gauge-dependent quantities — such as gluon and 
quark propagators — are evaluated, especially if 
the effect of Gribov copies is considered. 

The main issue regarding the efficiency of these 
algorithms is the problem of critical slowing-down 
(CSD), which occurs when the relaxation time 
t of an algorithm diverges as the lattice volume 
is increased (see for example ||). Conventional 
local algorithms have dynamic critical exponent 
z ~ 2, namely r ~ TV 2 , while global methods may 
succeed in eliminating CSD completely, i.e. z ~ 0. 
We consider five different algorithms: the Los 
Alamos method (a conventional local algorithm), 
the Cornell method, which is generally believed 
to have z w 2, the overrelaxation and stochas- 
tic overrelaxation methods (improved local algo- 
rithms, which are expected to show z » 1), and 
the Fourier acceleration method, which is a global 
method. We confirm the predictions for z (in 
both two and four dimensions) with the excep- 
tion of the Cornell method. For this case we ob- 
tain that the dynamic critical exponent is actu- 
ally z < 1, a result that can be understood by 
a comparative analysis between the Cornell and 
the overrelaxation methods. Besides the prob- 
lem of CSD, we are also interested in understand- 
ing which quantities should be used to test the 
convergence of the gauge fixing, and in finding 
prescriptions for the tuning of parameters, when 



needed. 

We consider the Standard Wilson action in d 
dimensions: 
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In order to fix the Landau gauge, we look for a 
local minimum of the function [Q 
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where 

uW{x) = g(x) U^x) g\x + ae M ) 



(2) 
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and we start from g(x) = I_ for all x, keeping the 
thermalized configuration {U^(x)} fixed. [For the 
SU (2) matrices we employ the usual parametriza- 
tion U = uo 1L + i u ■ a, where the components of 
a are the three Pauli matrices.] 

If {ujf (x)} is a stationary point of 8 ({<?}) 
then we have the well-known result 
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namely the lattice divergence of each color com- 
ponent of the gauge field 



A^x) = 
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is null. From equation (|J) it follows that 

q v {x v ) = j2J2M x ) 

is constant, namely it is independent of x„. 
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The algorithms we consider are all (with the 
exception of the Fourier acceleration) based on 
updating a single-site variable g(y) at a time. As 
a function of g(y) only, the minimizing function 
becomes 

a d 

£ \g(y)] = constant - tr[g(y) h{y)] , (7) 

where 

+ Ul(y - ae„) g\y - ae„) ] (8) 

is the single-site effective magnetic field. Note 
that h(y) is proportional to an SU(2) matrix, i.e. 
we can write 



h(y) = M{y)h{y) = y/ det h(y) h(y) 



(9) 



We also define T(y) = tr[g(y) h(y)]. 

The single-site (multiplicative) update can be 
written as 
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and for the methods we consider we have (sec |ij 
and references therein): 

Los Alamos Method: in this case we have 
g( new \y) — fv(y), i.e. this update brings the 
single-site function £ [g(y)] to its unique absolute 
minimum; 

Overrelaxation: here the matrix ft( u P date ) (y) 
is given by [ht (y) (y)] u , with id G (1,2); notice 
that the Los Alamos method corresponds to the 
case in which to is equal to one, and that for u> = 2 
the value of £ [g{y)] does not change; 
Stochastic Overrelaxation: in this case the 
update g^ new Hy) is given by 

h\y)gHy)"\ g(y) = h(y)T(y)- g(y) 

with probability p, and by h'(y) with probability 
1-p; 

Cornell Method: here Ri u P date ) (y) is propor- 
tional to [L-aa 2 g (V • A^) (y)] ; in § we 
prove that this can also be written as 



aM{y)T(y) 



lL + aM(y)rf{y)gi(i 



Fourier Acceleration: in this case the matrix 
j^(update) ^ j g p r0 p 0r tional to 
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where F is the Fourier transform, p 2 is defined 
as (4/a 2 ) sin 2 (nak^), and ak^N takes the 
values 0, 1, ...,N — 1. (Of course, to reduce 
the number of times the Fourier transform is 
evaluated, a checkerboard update should be em- 
ployed.) 

In order to determine the convergence of the 
gauge fixing, we measure the relaxation times Tj 
for the following quantities (lj: 



ei(t) 
ea(*) 

e 4 (i) 
e 6 (t) 
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Here t indicates the number of sweeps of the lat- 
tice and Q v = 1/N J2x =i Qv{ x i>)- We expect 
to observe that 



ei (t) oc exp ( — t/rj ) with 



dN** 



(11) 



and also that all the quantities above have the 
same relaxation time r, and hence the same c, z. 

2. RESULTS FOR d = 2 

For the four local methods we used Q lattice 
sizes N = 8, 12, ... , 36, while in the Fourier ac- 
celeration case we considered N = 8, 16, 32, 64. 
In all cases we have chosen the constant physics 
N 2 /(3 = 32, namely N/£ w 7. We stopped the 
gauge fixing when the condition e2 < 10 -12 was 
satisfied. 

From our data it is clear that the four quanti- 
ties ei have the same r for each given algorithm, 
as expected. However, for all the local updates 
- with the exception of the stochastic overre- 
laxation — the quantity ee was hardest to re- 
lax, namely its value was several orders of magni- 
tude larger than the value of the other quantities. 
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On the contrary, the Fourier acceleration method 
seems to be very efficient in relaxing e@- Basing 
on these results, we think that the quantity eg 
allows a very sensible check of the gauge fixing 
and, in our opinion, it should always be used. 

We evaluate the dynamic critical exponents z 
from the weighted least-squares fit for r = cN z , 
using lattice sizes N > N min , and obtain |Q: 



Los Alamos 


z 


= 1.99 ±0.04 , 


min 


= 12 


Cornell 


z 


= 0.83 ±0.09 , 


N ■ 

1 Y mm 


= 16 


Overrelax. 


z 


= 1.12 ±0.07 , 


^* min 


= 16 


Stoch. Overr. 


z 


= 1.09 ±0.05 , 


■N' min 


= 12 


Fourier 


z 


= 0.04 ±0.06 , 


N> min 


= 8 



Clearly, the Fourier acceleration method is the 
most successful in reducing CSD, and it should 
be the method of choice when large lattice sizes 
are employed. However, the method is more 
costly per iteration than the local methods, and 
its CPU-time/site increases logarithmically with 
the lattice size. When this is taken into account, 
the total time for gauge fixing a configuration is 
smaller for Fourier acceleration than for the local 
methods if the lattice size is larger than about 
300 sites, as can be seen in Figure 1. [Of course 
this result is very machine- and code-dependent. 
In any case, it seems unlikely jl[ that the Fourier 
acceleration method would become the method of 
choice in two dimensions at lattice sizes N smaller 
than around 100 sites.] 

3. TUNING 

The need for tuning the parameters of an al- 
gorithm in order to get optimal efficiency is, of 
course, a potential disadvantage. Of the five al- 
gorithms we consider, all (but the Los Alamos 
method) require tuning. We did a careful analy- 
sis of this problem and we were able to verify a 
simple analytic expression for the optimal choice 
of u> (overrelaxation method), and to relate u op t 
to the optimal choice for the parameters of the 
other local methods. More specifically: 
Overrelaxation: We tested the Ansatz || 

u opt = 2 [1 + Copt/NT 1 (12) 

and obtained very good agreement for lattice sizes 
greater than 12 with C op t — 1.53 ± 0.35. Note 



T„j for Cornell and Fourier acceleration 




Figure 1. Comparison of the "effective" CPU- 
time/site (T g f) between the Cornell method (the 
nearly straight line) and the Fourier acceleration. 



that, as the number of iterations t increases, 
the matrix R( u P date ) should approach the iden- 
tity matrix 1L. In this limit, the overrelaxation 
update can be written as jlj 

g {over \y) oc (l-w)ff(v) + w#(y). (13) 

Cornell Method: In this case, in the limit of 
large t, we can write Q 

g {Corn \y) oc [l-aM{y)]g{y) + aN{y)~h\y) . 

By comparing this equation with (|l3|), we made 
the conjecture a opt (A/ - ) = uJ op t, with (A/") given 
by 2d (I — {Smin}) ■ This Ansatz is very well 
fitted by our data 

Stochastic Overrelaxation: In this case it is 
not clear how to write a formula for the update as 
the number of iterations t increases. A possibility 
|Q] is to consider 

g^ oc \y) « (l + p)hHv)-pg(v), (14) 

which suggests the relation p op t — uj op t — 1. This 
conjecture seems to be satisfied reasonably well 
for lattice sizes larger than 30. 
Fourier Acceleration: By analogy with the 
continuum U(V) case, we expect [|| 

a 2 g (V-A^)(y)oiexp[-ata 2 pl ax ] , (15) 
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and therefore with the choice a op t = (ap max )~ 2 
we should obtain r oc 0(1). This gives, in d = 
2, a value a op t = 0.125. This result is only in 
qualitative agreement with our data, which seem 
to indicate that a op t ~ 0.16 for large lattice sizes. 

4. EXTENSION TO d = 4 

We study (in two and in four dimensions) the 
case j3 — oo, namely we fix U^(x) = 1L in (||) 
for all x and /i, and we look for a local minimum 
of the function (^]) starting from a random con- 
figuration {g(x)}. In both the two and four di- 
mensional cases, we essentially confirm the results 
obtained for z at finite (3 and in two dimensions 

It is interesting to notice that, in the case j3 — 
co, we can also write the minimizing function as 



£({<?}) 



2dV ^ 

x, n 



\g(x) - g{x + ae„) f. (16) 



Here the SU(2) matrix g = g^lL + i 3 ■ g is con- 
sidered as a four-dimensional unit vector (go,g) 
and ||g|| 2 = g\ + g ■ g. When we are close to a 
minimum we can write 



g(x) = ±-iea-f{x) + 0(e 2 ) , 



(17) 



where f(x) is a three- vector field. This gives 
g{x)g\x) = i_ + 0(e 2 ) (18) 



and 



£({g}) 



2dV 



^||/(x) - /(a; + aeAl )|| 2 ,(19) 



x, ll 



namely (up to order e 2 ) we have the action of 
a three- vector massless free field f{x). Therefore 
we can use standard analytic methods in order 
to study the problem of minimizing this quadratic 
form, and to compare the results for z, t and the 
optimal choice of the parameters with our numer- 
ical data. Indeed, we find good agreement for all 
the quantities and methods considered j| . 

Finally, we extend our simulations to the case 
(3 = (again in two and in four dimensions) gj. 
All our results at (3 > are confirmed except 
for the Cornell method, for which we get z ss 2, 



and for the Fourier method, which shows z w 
1. In the Cornell case this result can be easily 
understood, since at (3 = its single-site update 
does not decrease the minimizing function at each 
step. In particular, the value of £ increases if 
aJV(y) is large than two. [This is obvious if we 
consider the analogy between the Cornell method 
and the overrelaxation method, and the relation 
aN(y) ~ ui. Moreover, it is plausible that only at 
very small values of (3 we can have a (Af) smaller 
than two and, at the same time, aj\f(y) > 2 for 
a large set of lattice sites y] However, we can 
recover the value z « 1 by redefining the Cornell 
single-site update in the following way || 
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(Corn) 



(y) OC 



1 



a(y)T(y) 



g{y)+a{y) h){y) , 



with a(y) = min ( aj\f(y), 2 ). In the Fourier ac- 
celeration case we have the same problem: the 
update can increase the value of the minimizing 
function £, and this happens most likely for small 
values of (3. However, in this case we did not 
find a redefinition of the Fourier update H which 
could recover the value z ~ 0. 
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